APS/123-QED 



Susceptibilities and fluctuations in a Quark-Hadron System 
with Dynamical Recombination 

Stephane Haussler and Marcus Bleicher 
Institut fur Theoretische Physik, Johann Wolfgang Goethe-Universttdt, 
Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany 



oo 

o 

O 

(N 



in: 

CD 



> 

00 

cn 

o 

00 
O 



X 



Horst Stocker 

Gesellschaft fur Schwerionenforschung ( GSI), Planck-Str. 1, D-64291 Darmstadt, Germany 
Institut fur Theoretische Physik (ITP), Johann Wolfgang Goethe- Universitat, 
Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany and 
Frankfurt Institute for Advanced Studies (FIAS), 
Ruth-Moufang-Str. 1, D-60438 Frankfurt am Main, Germany 

(Dated: March 19, 2008) 

Within the framework of the dynamical recombination approach implemented in the par- 
tonic/hadronic quark Molecular Dynamics (qMD) model, we investigate numerous QGP signals 
constructed from the correlations and fluctuations of conserved charges, namely charged particle 
ratio fluctuations, charge transfer fluctuations as well as various ratios of susceptibilities. We argue 
that more generally, the covariances and the variances of the conserved charges divided by the num- 
ber of charged particles are a measure of the quark number susceptibilities and are thus sensitive 
to the phase of the system. Computations carried within samples of central qMD events from low 
AGS energies on (Ei a t = 2 AGeV) up to the highest RHIC energies available (^/sjvjv = 200 GeV) 
show that the final state calculations are always compatible with the hadronic result. From anal- 
yses performed as a function of time with our sample of events at ^snn = 200 GeV, we find 
the recombination-like hadronization process the qMD system undergoes to be responsible for the 
vanishing of the early stage QGP correlations and fluctuations. These results might explain the 
compatibility of the measurements on charged particle ratio fluctuations with the hadronic expecta- 
tions and leave no room for the experimental observation of any QGP fluctuation/correlation signals 
based upon the susceptibilities of the conserved charges. 

PACS numbers: 25.75.Nq,24.60.-k,12.38.Mh 



I. INTRODUCTION 

There exist a large consensus in the heavy-ion physics 
community that a plasma of deconfined quarks and glu- 
ons is formed during the first instant of the collision 
of heavy nuclei for sufficiently high incoming energies 
[E 0- H| ■ The investigation of the properties of the highly 
excited nuclear matter produced requires the use of phys- 
ical quantities sensitive to the initial partonic state of the 
system which survive hadronization and the subsequent 
hadronic evolution. 

Numerous probes based on the eventwise fluctuations 
of physical quantities have been proposed to_inve_stigate 
the state and the properties of the fireball [ 
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them, charged particle ratio fluctuations, charge transfer 
fluctuations, baryon number-strangeness correlations and 
charge-strangeness correlations were prominently pro- 
posed to pin down the formation of the deconfined phase 

at rhic (HI, H2, HE HE H3, m, m, SO, [U, [43, El] . 

It was pointed out that these quantities should reflect the 
properties of the system in the first instants of the reac- 
tion and survive the whole evolution: with a strong lon- 
gitudinal flow, any conserved charge is frozen in a given 
rapidity window because the expansion is too fast for the 



initial partonic fluctuations to relax through the trans- 
port of charges in and out of a given rapidity slice. Thus, 
QGP fluctuations should not have the possibility to relax 
to their hadronic expectation values. It is clear that the 
size of the rapidity window for fluctuation studies must 
not be too wide in order to avoid global charge conser- 
vation which would lead to a vanishing signal, but also 
neither too small to avoid purely statistical fluctuations 
and the transport of charges in and out of the slice under 
consideration by hadronic rescattering. The generally ac- 
cepted rapidity width is of the order of Ay = 0.5 — 1 units 
in rapidity. 

However, hadronization itself might bring the initial 
QGP correlations and fluctuations to values compatible 
with the hadronic expectations. How to calculate or de- 
scribe the non-perturbative domain of QCD where the 
physics of hadronization sits remains today an open ques- 
tion. A possible mechanism consists of the recombination 
of quarks into hadrons 0, HE HE S3, EH, EE HE HI HI 
HE Hi]]. On the experimental side, the most striking evi- 
dences for recombination are the number-of-constituent- 
quarks scaling of the elliptic flow vi [IE HE HH and the 
large number of baryons over number of mesons ratio in 
the intermediate pt range measured in central Au+Au 
collisions at RHIC energies [IE HE Hz|- In this pa- 
per, we use the quark Molecular Dynamics (qMD) model 
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(HH HH, [H, H3| where quarks coalesce locally in coordi- 
nate and momentum space due to a confining inter-quark 
potential. 

In previous studies [H, [33[ performed with qMD and 
in an earlier exploratory work (58j |. it was argued that 
recombination leads to a complete vanishing of various 
initial state QGP fluctuation and correlation signals. We 
extend this idea further and show that dynamical recom- 
bination generally blurs QGP signals constructed from 
the correlations and fluctuations of conserved charges. 

Besides the predicted impossibility to experimentally 
establish the existence of QGP fluctuations, our results 
are important for the fluctuation investigations planned 
for the search of the critical endpoint of the QCD phase 
diagram at the forthcoming FAIR facility at GSI. 

The article is organized as follow: We first introduce 
the quark Molecular Dynamics (qMD) model and ad- 
dress the problem of entropy conservation at the par- 
tonic/hadronic transition which might occur in some re- 
combination approaches. We detail the situation for 
our dynamical recombination procedure and argue that 
we respect the necessary increase of entropy with time. 
Within our sample of events generated with the qMD 
model, we then investigate the charged particle ratio 
fluctuations, the charge transfer fluctuations, the baryon 
number-strangeness correlation coefficient, the charge- 
strangeness correlation coefficient and discuss the influ- 
ence of recombination-hadronization for all these observ- 
ables. We also compare the result of the qMD with lat- 
tice data [5^ . [60j on various ratios of susceptibilities and 
to directly investigate the susceptibilities by performing 
event-by-event analyses. Finally, we discuss the reasons 
for the transition from the initial QGP fluctuations and 
correlations to their hadronic value. 



II. THE qMD MODEL 

The qMD model 0, HI HE HI employed here 
is a semi-classical molecular dynamics approach where 
quarks are treated as point-like particles carrying color 
charges and interact via a linear heavy quark potential. 
Initial conditions for the qMD are taken from the hadron- 
string transport model UrQMD [f5l], [HJ : After the two 
incoming nuclei have passed through each other, (pre-) 
hadrons from the strings and fully formed hadrons from 
the UrQMD model are decomposed into quarks with cur- 
rent masses m u — md — 5 MeV and m s — 150 MeV. At 
the highest RHIC energy, this happens at a center of 
mass time of t = 0.15 fm/c. It should be noted that the 
qualitative results of the present study are not restricted 
to any specific initial state. The UrQMD model is solely 
used to provide an exemplary initial state after the ini- 
tial qq production has taken place. The quarks are then 
let to evolve and interact within the qMD via a linear 
potential Vflr^ — r,-|) = «|rj — r,-|, where k is the string 
tension and r„ is the position of particle n. Therefore 



the full Hamiltonian of the model reads: 
N 1 1 

« = E\/^ m . 2 + 2^ c ' 3y(|r, " rjl) ■ (1) 

where N counts the number of particles in the system and 
the term Cy takes into account the color dependence of 
the interaction. 

The quark-(anti-)quark interaction within this poten- 
tial naturally leads to confinement through the binding 
of (anti-)quarks into color neutral clusters. New hadrons 
are formed from quarks whose momentum and positions 
are close to each others. Typical values for the relative 
momenta of the quarks in the two-particle rest frame at 
hadronization are \p q \ = \pq\ < 500 MeV, the typical 
distance is below 1 fm, i.e. hadronization occurs locally 
into hadronic clusters of mesonic and baryonic type that 
resemble the Yo-Yo states of the LUND model. Most 
of the mass of the newly formed hadrons come from the 
potential energy between the quarks forming the con- 
sidered cluster. We generally obtain clusters with large 
masses mapped to hadronic resonances and thus describe 
the transition from current quark masses to constituent 
quark masses of the order of 300 MeV after hadroniza- 
tion has occured. These clusters are then mapped to 
hadron states that are later allowed to decay in the fur- 
ther evolution of the system. The reader is referred to 
[5TI . I52I l53l Hi[ for a detailed discussion of the qMD 
model. 



III. ENTROPY CONSERVATION 

The question of entropy at the parton/hadron tran- 
sition needs to be examined when quark recombination 
is used to describe hadronization. Entropy is propor- 
tional to the number of particles in a gas of massless par- 
tons/hadrons. Thus, when partons recombine to form 
new hadrons, the total number of particles decreases and 
with it entropy. This situation is not acceptable as it 
violates the second law of thermodynamics. 

Contrary to many other recombination models, one can 
argue that qMD does not violate the entropy condition. 
In fact, entropy does not only depend on the number of 
particles but also on their masses as examined in [58| . 
When m/T > 3, entropy can be simply approximated 
for a gas of massive particles through: 

S = Y(3.5 + m/T) , (2) 

where S stands for entropy, N for the number of particles 
in the system, m for their average mass and T for the 
temperature at the transition. 

Fig. Q] shows the distribution of the clusters masses at 
hadronization obtained from qMD central Au+Au sim- 
ulations at y/s — 200 AGeV. The average mass after 
recombination of the quarks is large with M c i uster ~ 
730 MeV. 
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FIG. 1: Distribution of the mass of the clusters obtained after 
recombination of the initial partons for central Au+Au qMD 
reactions at ^/snn = 200 GeV. 




FIG. 2: Total number of particles N qua rks + N^adrons (full 
circles) as a function of time and average mass as a function of 
time (open circles) for Au+Au reactions at ^/snn = 200 GeV. 



To estimate entropy, let us consider a system of N q+q 
quarks and anti-quarks that recombine into N q+q /2 
hadrons where baryons are neglected for simplicity. The 
quarks being nearly massless, entropy in the partonic 
phase is equal to [581 ] : 

S QGP = A.2N q+q , (3) 

After recombination and using Eq. [5] to estimate entropy 
in the hadronic phase: 

Sreco = * (3 - 5 + 730 /150) = 4.2N q+q . (4) 

where the qMD average cluster mass m = 730 MeV and 
the qMD critical temperature T = 150 MeV are used. 
Thus, SVeco ~ Sqgp and entropy does not decrease in 
the process of recombining quarks into massive hadronic 
clusters. 

Entropy is further analyzed through the number of par- 
ticles in the qMD model. Fig. [2] depicts the total number 
of particle (full circles) and the average particle mass 
(empty circles) as a function of time. The average parti- 
cle mass peaks around a time t s» 11 fm/c to an average 
mass of (to) ~ 450 MeV. This is much less than the av- 
erage cluster mass observed in Fig. [1] and is due to the 
decay of the excited clusters produced at hadronization. 
Further decrease of the average particle mass is related 
to the continuing decay of resonances when all particles 
are already confined into hadronic clusters. As can be 
seen from Fig. [21 the total number of particles decreases 
at hadronization. However, resonance decay leads to a 
number of particles in the final state close to the initial 
number of quarks. The increase of entropy with time is 
then not violated, at least when comparing initial and 
final state. 



Results from lattice QCD calculations [H| also indicate 
that recombination is compatible with entropy conserva- 
tion at the transition. 



IV. FLUCTUATIONS OF CONSERVED 
CHARGES 

Let us set the stage by exploring the time evolution 
of the hadronization dynamics in the model. Fig. de- 
picts the fraction of quark matter on the total num- 
ber of particles in the system, i.e. quark fraction = 
(n q + n q ) I (rihadron + n q + n q ) as a function of time for 
central Au+Au qMD events at ^/sWn = 200 GeV. With 
the given initial conditions, the fireball stays in a decon- 
fined state during the first 6 fm/c where almost no quarks 
hadronize. As the system expands further and density 
decreases, quark recombination into baryons and mesons 
occurs and the number of deconfined quarks drops to 
zero. The point at which the parton-hadron transition 
occurs will be denoted by arrows in the following plots. 
We now turn to the investigation of the various fluctua- 
tion signals. 



A. Charged Particle Ratio Fluctuations 

Charged particle ratio fluctuations were proposed as 
a clear signal for the onset of the quark-gluon plasma 
phase 39] . The basis for the argument is that the quanta 
of the electric charge carrier are smaller in a quark gluon 
plasma phase than in a hadron gas and are distributed 
over a larger number of particles. Moving one charged 
particle from/to the rapidity window then leads to larger 
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FIG. 3: Fraction of the number of quarks from the total num- 
ber of particles (N q + Nq)/(N q + Nq + Nh) as a function 
of time at midrapidity for central Au+Au qMD reactions at 
^/sjvjv = 200 GeV. The arrow depicts the hadronization time. 

fluctuations in a hadron gas than in a QGP. The charged 
particle ratio fluctuations can be quantified by the mea- 
sure D calculated in the rapidity window Ay: 

D = -^ 7r (N ch )(6R 2 ) Ay . (5) 

where N c h stands for the number of charged particles, 
R = (1 + F)/(l - F) with F = Q/N ch , Q being the 
electric charge. Following [l5|, the uncorrected charged 
particle ratio fluctuations D = (N c i l )(8R 2 ) are divided 
by the factors C M to correct for finite net charge due to 
baryon stopping and C y for global charge conservation. 
It was argued that depending on the nature of the initial 
system, D will yield distinctly different results: D = 1 
for a quark-gluon plasma, D — 2.8 for a resonance gas 
and D = 4 for an uncorrelated pion gas. 

Experimentally, charge fluctuations have been mea- 
sured at RHIC energies by STAR [H, |g| and PHENIX 
[6?1 [67J ■ Both experimental analyses yield results com- 
patible with a hadron gas. Further results from CERN- 
SPS [H,[6^| based on a slightly different measure are also 
compatible with the hadronic expectations. 

Fig. 2] depicts the energy dependence of the corrected 
charged particle ratio fluctuations D from AGS ener- 
gies on (Ei a t,=2 AGeV) up to the highest RHIC en- 
ergy available (y/~SNN=200 GeV) calculated with central 
Au+Au/Pb+Pb qMD events. As suggested in [IE El, 
D is calculated in a rapidity window of y — ±0.5. D de- 
creases from D » 6 down to D 3 with increasing energy 
and is not compatible with the QGP value in the whole 
energy range explored, even though the system initially 
went through a quark phase. No hadronic rescattering 
stage being included in the present calculations, this is a 
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FIG. 4: Corrected charged particle ratio fluctuations D as a 
function of the center of mass energy ^/sjvjv within the qMD 
model for central Au+Au/Pb+Pb reactions (full symbols). 
Also shown are the values for an uncorrelated pion gas and a 
quark-gluon plasma. 



first indication that the hadronization procedure imple- 
mented in the qMD model alone destroys the early stage 
QGP fluctuations. 

The signal might also be sensitive to the size of the 
rapidity window A y u sed as can be seen from VNI 
model calculations [41] . The rapidity window depen- 
dence for both corrected D (full symbols) and uncor- 
rected D (empty symbols) charged particle ratio fluctua- 
tions is plotted on figure Fig. [5] for central Au+Au qMD 
events at ^/snn — 200 GeV where the quark phase du- 
ration is the longest in the present calculations. One can 
measure the importance of the correction terms as the 
rapidity window is widened and goes to a regime where 
global charge conservation becomes so important that D 
drops to nearly D 0. When the rapidity window is 
wide enough to take into account all particles present in 
the events, the net charge is fixed by the initial charge of 
the colliding gold nuclei and can not fluctuate any more. 
Furthermore, both D and D come very close to 4 when 
the rapidity window Ay is decreased to very small values. 
If the rapidity window is so small that only one of the 
decay products of a resonance can be observed, then the 
system appears completely uncorrelated. In the present 
model, one observes that D decreases from D = 4 down 
to D w 3 when Ay w 1 and then does not change with 
increasing rapidity window width until Ay « 3. The re- 
sult always stays at the expected hadronic expectation 
value. Thus, the qMD result does not depend on the size 
of the rapidity window used and is compatible with the 
hadronic expectations. 

Within the simulations, one has access to the evolution 
of D from the initial quark phase to the final state. Fig.[S] 



5 



5 
4.5 



Uncorrelated Pion Gas 




0.5 



1.5 

Ay 



2.5 




10 

time [fm/c] 



FIG. 5: Corrected D (full symbols) and uncorrected D (empty 
symbols) charged particle ratio fluctuations as a function of 
the width of the rapidity window within the qMD model for 
central Au+Au reactions at -^/sjvjv = 200 GeV. Also shown 
are the values for an uncorrelated pion gas, a resonance gas 
and a quark-gluon plasma. 



FIG. 6: Corrected charged particle ratio fluctuations D = 
D/(CyCft) as a function of time within the qMD model for 
central Au+Au reactions at y'sjvjv = 200 GeV (full symbols). 
Also shown are the values for an uncorrelated pion gas, a 
resonance gas and a quark-gluon plasma. The arrow depicts 
the hadronization time. 



shows the result for D from the qMD recombination ap- 
proach as a function of time calculated for a sample of 
central Au+Au events at ^snn = 200 GeV in a rapidity 
window of y = ±0.5. In the early stage, when the system 
is completely in the deconfined phase, D as 1 as expected. 
When approaching the hadronization time, D starts to 
increase and reaches D as 3.5 after hadronization. The 
increase of D occurs exactly at the same time as the re- 
combination of the quarks and anti-quarks into hadrons. 
The slight decrease of D at later times is related to the 
continuing decay of resonances lowering the charge ratio 
fluctuations @. 

Gluons are not explicitly included in our simulations 
and their influence at the transition should be discussed. 
Following [39(, the uncorrected charged particle ratio 
fluctuations D can be expressed as: 



D 



(SQ 2 ) 
' (Nch) 



(6) 



where Q stands for the net charge and (N c h) for the mean 
number of charged particles. Gluons neither participate 
in the fluctuations of the charge nor to the number of 
charged particles in the pure deconhned phase. However, 
when the system hadronizes, gluons participate in the 
denominator of Eq. [5] through gluon fragmentation: g — > 
hadrons, which can lowers D. In qMD, the decay of 
highly excited hadronic clusters, whose masses include 
the inter-quark potential, effectively models a part of this 
contribution. 



B. Charge Transfer Fluctuations 

As a next observable, we turn to charge transfer fluctu- 
ations that were also suggested to provide insight about 
the formation of a QGP phase. Charge transfer fluctua- 
tions are a measure of the local charge correlation length. 
They are defined as [H, El, [H] : 



D U ( V ) = (u(r,) 2 



(«fa)> S 



(7) 



with the charge transfer u(rj) being the forward-backward 
charge difference: 



u(r)) = [Q F (r,)-Q B ( V )]/2 , 



(8) 



where Qf and Qb are the net charges in the forward 
and backward hemisphere of the region separated at n 
and calculated in the pseudo-rapidity window A?/ = r/± 
1. 

The charge correlation length and hence the local 
charge fluctuations are expected to be much lower in a 
quark-gluon plasma than in a hadron gas. Because the 
measured quantity is local, it can provide information 
about the presence and the extent of a QGP in (pseudo- 
rapidity space. Thus, one expects to observe the lowest 
value of the charge transfer fluctuations at midrapidity, 
where the energy density is the highest and where the 
plasma should be located. Experimental data on this 
observable are not available up to now. 

The time evolution of the charge transfer fluctuations 
from the present calculations are shown in Fig. [7] The 
sample of events consists again of central Au+Au reac- 
tions simulated at ^/snn = 200 GeV. Here, the pseudo- 
rapidity window corresponds to the STAR acceptance 
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FIG. 7: Charge transfer fluctuations at midrapidity (A77 = 
±1) and separation at 77 — as a function of time within 
the qMD model (full symbols) for central Au+Au reactions 
at t/snn = 200 GeV. The arrow indicates the hadronization 
time. 



FIG. 8: Charge transfer fluctuations for a total rapidity win- 
dow Ar/ = 77 ± 1 as a function of pseudo-rapidity r\ with 
the qMD model for central Au+Au reactions at ^/sjvjv = 
200 GeV. 



(77 = ±1). As expected, the correlation length (at central 
rapidities) is small, with D u j \dN c h/dy) » 0.1, as long 
as the system is in the quark phase. However, similar to 
the charged particle ratio fluctuations discussed above, 
the charge transfer measure increases with time up to its 
hadronic value of D u /(dN c h/dy) w 0.5 when hadroniza- 
tion has occured. The final state result is in agreement 
with the value given by HIJING calculations and there- 
fore in line with the hadronic expectation [3^, HU, Hf| HI] • 

The dependence of charge transfer fluctuations 
D u / (dN c hldy) on pseudo-rapidity calculated from a sam- 
ple of central Au+Au qMD events at ^Jsnn = 200 GeV 
is depicted on Fig. [51 The pseudo-rapidity dependence of 
this observable is flat and compatible with the hadronic 
expectations (D u /(dN c h/dy) w 0.5). In particular, no 
dip is visible at 77 ~ where the plasma is located. There- 
fore, charge transfer fluctuations seem not only insensi- 
tive to the initial QGP distribution in pseudo-rapidity 
space, but neither to the existence of the quark phase 
itself in our simulations. 

Fig. [HI depicts the excitation function of D u /(dN c h/dy) 
from central Au+Au/Pb+Pb qMD events calculated at 
midrapidity in the pseudo-rapidity window 77 = ±1. 
D u / (dN c h/dy) increases from D u /(dN c h/dy) » 0.3 at 
low AGS energies to reach D u /{dN c h/dy) s» 0.5 at RHIC 
energies. The result stays above the QGP expectations 
in the whole energy range explored. The qMD estimates 
for charge transfer fluctuations are consistent with the 
model studied in [70l |. also including a partonic/hadronic 
transition. 
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FIG. 9: Charge transfer fluctuations at midrapidity (A77 = 
±1) and separation at 77 = as a function of the center of 
mass energy y / ijvjv within the qMD model (full symbols) for 
central Au+Au/Pb+Pb reactions 



C. Ratios of Susceptibilities 

We now turn to the study of the baryon number- 
strangeness correlation coefficient Cbs, charge- 
strangeness correlation coefficient Cqs and more 
generally to the ratios of susceptibilities. 

It was argued in 38] that the correlation between 
baryon number and strangeness strongly differs whether 
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the system is in a QGP or in a hadronic phase. In 
a QGP, strange quarks will always carry baryon num- 
ber ±i, which strongly correlates baryon number and 
strangeness. The situation is different in a hadron gas 
where strangeness can be carried without baryon num- 
ber, e.g. with kaons. Following [38|, we calculate the 
Cbs correlation coefficient on an event-by-event basis: 



Cbs 



, (BS)-(B)(S) 



(9) 



where B stands for the baryon number and S for the to- 
tal strangeness. The expectation value in an ideal QGP 
is Cbs — 1 and Cbs = 0.66 in a resonance gas at tem- 
perature T = 170 MeV and chemical potential /i = 
corresponding to the conditions obtained at RHIC 1 38] . 

Due to difficulties concerning the measure of neutral 
particles, Cbs is difficult to extract from the experimen- 
tal data. It is however possible to use the related cor- 
relation coefficient Cqs, where Q stands for the electric 
charge (59|: 



Cqs 



, (QS)-(Q)(S) 



(10) 



Under the assumption of isospin symmetry [60| . Cbs 
and Cqs are related through: 



Cqs — 



3 — Cbs 



(11) 



The variances and covariances of conserved charges 
can also be expressed in term of susceptibilities and off- 
diagonal susceptibilities x- 



((8Xf) = VT XX , 
((SX)(6Y)) = VT X xy , 



(12) 



where V is the volume and T the temperature. 

Taking the ratio of different susceptibilities then makes 
the volume term disappear. The Cbs and Cqs coeffi- 
cients can thus be written as: 



Cbs = 



Xs 



Xs 



(13) 



Fig. [10] (left) depicts Cbs and Cqs as a function of 
the temperature as measured by the kinetic energy of 
the partons in the qMD model (full symbols). In com- 
parison, lattice QCD calculations ^59j (open symbols) are 
shown. At high temperature, 1QCD yield the results of 
an ideally weakly coupled QGP with Cbs = Cqs — 1 
while below T c , they are consistent within uncertainties 
with a hadron gas. qMD calculations yield results with 
Cbs ~ 0.7 and Cqs ~ 1-2 in the low temperature regime. 
Both quantities also go to Cbs ~ Cqs ~ 1 in the quark 
gluon plasma phase. Similarly to charged particle ratio 
fluctuations and charge transfer fluctuations, the QGP 
signal vanishes in the hadronic phase. 

Various other ratios of susceptibilities using the electric 
charge Q, strangeness £>, baryon number B, third com- 
ponent of the isospin 1% and hypercharge Y have been 



measured on the lattice and can be used to distinguish 
between different models, or might serve as signals of the 
deconfined phase. 

Here, we show the result of xb/xq, Xh/XQ, Xy/xq 
and — (xs/xq)/2 as a function of temperature from 
qMD calculations (full symbols) and compare with the 
available lattice data (open symbols) on Fig. [10] (right). 
Xh/XQ is essentially flat with XI3/XQ ~ 0-9 with no dif- 
ference between the deconfined and the confined phase. 
On the contrary, all other quantities increase from the 
confined to the deconfined phase with xy Ixq going from 
Xy/xq ~ 0-6 to xy/xq ~ 0.7 in the temperature range 
explored. Xb/xq and — (xs/xq)/2 are roughly the same 
and increase from xb/xq ~ -(xs/xq)/ 2 ~ 0.2 to 
Xb/xq ~ -(Xs/Xq)/ 2 ~ 0.4 with increasing temper- 
ature. 

Even though the qMD approach is a grossly simplified 
version of QCD, it compares well with the lattice predic- 
tions over the whole temperature range accessible in the 
course of a collision. This might be seen as an a poste- 
riori justification for the use of the qMD model for the 
fluctuation and correlation analyses performed above. 



V. SUSCEPTIBILITIES 

Dividing the variances and the covariances of the con- 
served charges by the average number of charged parti- 
cles (N c h) (see Eq. IT2"j) . which is proportional to UT 3 , 
it is possible to access directly information on the diag- 
onal and off-diagonal susceptibilities x- We define the 
coefficient oxy by : 



&XY 



(XY) 



(X)(Y) xxy 
I) 



(14) 



where (X,Y) = {B,Q,S,I 3 ,Y, . . .). 

Note that ctqq corresponds to the net charge fluctua- 
tions and is equivalent to the charged particle ratio fluc- 
tuations investigated before (see Eq. [S]). 

Fig. [TTJ (left) depicts the time dependence of oqq (full 
squares), ass (open squares) and obb (full circles). All 
these quantities exhibit a sharp rise in the vicinity of 
the hadronization time followed by a smooth decrease to 
their final hadronic value: (Jqq increases from uqq ~ 0.2 
to tqq ~ 0.7, consistently with the results obtained for 
charged particle ratio fluctuations. <jbb increases from 
obb ~ 0.08 to obb ~ 0.12 and ass from obb ~ 0.14 to 
(J bb ~ 0.24. Besides the notable exception of ctqq, they 
all exhibit a characteristic enhancement at the transition 
from partonic to hadronic matter. 

The result of the correlation coefficients <jxy is pre- 
sented in Fig. [TTJ (right) for uqs (full squares), ass (open 
squares) and <jqb (full circles), ass is multiplied by the 
factor -1. The behaviour is qualitatively similar with uqb 
increasing from oqb ~ 0.02 in the pure deconfined phase 
up to oqb ~ 0.03 in the hadronic stage. <jqs increases 
with time from <jqs ~ 0.05 up to <jqs ~ 0.1. <jbs ~ 0.05 
in both partonic and hadronic phase but still exhibit an 
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FIG. 10: Left: Cbs and Cqs correlation coefficient as a function of temperature calculated with the qMD model (full symbols) 
and comparison to lattice data (empty symbols). Right: Various ratios of susceptibilities Xx/xq except for X = S which 
is (xs/xc?)/2- Open symbols are the result of lattice calculations [591 ]. Full symbols are the result of the qMD model. The 
temperature is extracted out of the kinetic energy per quark in the medium. 



0.8 
0.7 
0.6 
0.5 



c X0.4h 



0.3 
0.2 
0.1 





10 

time [fm] 



10 2 



0.14 
0.12 
0.1 

0.08 

>- 

x 

) 

0.06 
0.04 
0.02 





10 

time [fm] 



10 2 



FIG. 11: Left: ctqq (full squares), ass (open squares) and obb (full circles) as a function of time. Right: oqs (full squares), 
o"s5 (open squares) and oqb (full circles) as a function of time, ass is multiplied by the factor -1. Calculations are performed 
with a set of central Au+Au qMD events at ^/s = 200 AGeV. 



enhancement around the transition time. The situation 
for axY strongly contrasts with the ratios of susceptibil- 
ities which smoothly increase from their QGP to their 
hadronic values. 

These correlation and fluctuation studies of uxy can 
be extended to isospin 1% and hypercharge Y from 
which experimentally measurable quantities can be con- 
structed, e.g. the Cms correlation coefficient [6(3]. How- 
ever, the investigation of B, Q and S is sufficient to treat 
the problem completely for three light flavors of quarks 



u ,d and s. In the QGP stage, Q, B and S can be rewrit- 
ten in terms of net upness Au, net downness Ad and net 
strangeness As [5^| whose fluctuations and correlations 
translate into quark susceptibilities. Note that equivalent 
sets of conserved charges can be chosen, e.g. (B, Is,Y) 
and extended to heavy quark flavours. 
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FIG. 12: Mean number of charged particles (N c h) as a func- 
tion of time. Calculations are performed with a set of central 
Au+Au qMD events at yfs = 200 AGeV with resonance decay 
included (full circles) and without (open circles). 



VI. VANISHING OF THE SIGNALS 

Finally, we discuss the transition from the initial QGP 
fluctuations and correlations to their hadronic values for 
the observables calculated here. We first investigate the 
effect of the decreasing mean number of charged par- 
ticles (N c h) at hadronization entering the definition of 
charged particle ratio fluctuations, charge transfer fluc- 
tuations and <txy coefficients. 

To pin down the effect of the dynamical recombina- 
tion process itself, free from the influence of resonance 
decay, we use a set of central Au+Au qMD events at 
y/SNN — 200 GeV where the hadronic clusters formed 
are not allowed to decay. The time dependence of (N c h) 
is depicted on Fig. [12] with resonance decay (full cir- 
cles) and without resonance decay (open circles). In 
the absence of decaying hadronic clusters, the num- 
ber of particles decreases from {N c h) ~ 11000 down to 
{Nch) ~ 2700, whereas it decreases from (N ch ) sa 11000 
down to {N c h) ~ 5500 for the default calculations. The 
initial N q charged quarks recombine into N q /2 hadrons, 
among which half, i.e. N q /4, are charged. The variances 
and covariances of the conserved charges exhibit no sharp 
structure at the partonic/hadronic transition. Thus, the 
decrease of the number of charged particles at hadroniza- 
tion is responsible for the sharp jump of the oxy coeffi- 
cients, charged particle ratio fluctuations D and charge 
transfer fluctuations D u / (dN c h / dy) . For full qMD calcu- 
lations, resonance decay affects both the numerators and 
denominator of the studied observables and slightly low- 
ers down the result to values which are still compatible 
with the hadronic value. 

For the ratios of susceptibilities such as Cbs, this ar- 
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FIG. 13: Distribution of the rapidity shift of the quarks at 
hadronization y quar k — Y c i uster - The mean rapidity shift is 

(\yquark VeiusteT*!/ — 0.57. 



gument is not valid as {N c h) plays no role. The tran- 
sition from the QGP to the hadronic values can result 
from correlations induced by the interaction potential it- 
self, i.e. the gathering of quarks into color neutral clus- 
ters, the shuffling of the quarks in and out of the consid- 
ered rapidity window at hadronization, or be affected by 
resonance decay. To distinguish between these different 
mechanisms, we use sets of central Au+Au qMD events 
at y/sxN — 200 GeV with no mapping to hadronic clus- 
ters -the quark system is let to evolve according to the 
qMD Hamiltonian of Eq. Q3- to set apart the effect of the 
interaction potential, with mapping and no resonance de- 
cay to investigate the rapidity shift of the quarks from 
their initial value to the hadronic clusters rapidity, and 
with full calculations including resonance decay. 

Let us discuss the capability of the qMD model to 
move charges in and out of the studied rapidity win- 
dow at hadronization. Fig [T2] depicts the distribu- 
tion of the rapidity shift of the quarks at hadronization 
Vquark — ^cluster obtained in the qMD model. The mean 
rapidity shift is {\y quark -Y c i uste r\) = 0.57 and can partly 
explain the vanishing of the signals. 

Figs IT4l shows the time evolution of Xq/xs (top left), 
Xq/Xb (top right) and Xb/xs (bottom) from qMD cal- 
culations where the quark system is let to evolve accord- 
ing to Eq. [1] (full circles), when the quark cluster are 
mapped to hadronic states but resonances not allowed 
to decay (open squares) and finally for full qMD calcula- 
tions with mapping and resonance decay included (open 
circles) . 

The effect of the rapidity shift depicted in Fig. [13] is 
the difference between the calculations without mapping 
(full circles) and without resonance decay (open circles). 
Xq/Xs without mapping stays flat around xq/xs ~ 1-3 
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FIG. 14: Top Left: Xq/xs as a function of time Top right: Xq/Xb as a function of time. Bottom: xb/xs as a function of time. 
Calculations are performed with sets of central Au+Au qMD events at y^s = 200 AGeV without hadronization (full circles), 
with mapping and no resonance decay (open circles) and for full qMD calculations (open squares). 



whereas Xq/xs increases with time from xq/xs ~ 1-3 
U P to Xq/xs ~ 1-8- The further increase to xq/xs ~ 2.9 
in full qMD calculations is related to the decay of reso- 
nances, xq/xb increases smoothly from xq/xb ~ 2.3 up 
to xq/xb ~ 3 without mapping and from xq/xb ~ 2.3 
U P 1° Xq/xb ~ 3.8 with mapping. The further increase 
to Xq/xb ~ 6 is due to the continuing decay of reso- 
nances. The situation is markedly different for Xb/xs 
which decreases with time from Xb/xs ~ 0.57 down to 
Xb/xs ~ 0.47 with no visible difference between all sets 
of events. In this last case, the decrease of Xb/xs with 
time is related to the potential-driven gathering of the 
quarks into color neutral clusters. We thus conclude that 
the ratios of susceptibilities is influenced by these three 



different effects with varying importance, the correlation 
due to the inter-quark potential used, the rapidity shift 
of the quarks to the rapidity of the hadronic cluster at 
hadronization, and the decay of resonances. 



VII. CONCLUSIONS 

We have investigated various quantities constructed 
from the fluctuations and correlations of the conserved 
charges within the dynamical recombination model quark 
Molecular Dynamics (qMD) which includes an explicit 
microscopic transition from partonic to hadronic mat- 
ter. The sample of events used consists of cen- 
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tral Au+Au/Pb+Pb reactions simulated with the qMD 
model from low AGS energies on (Ei ab = 2 AGeV) up to 
the highest RHIC energy available (^/satjv = 200 GeV). 

We first discussed our approach and addressed 
the problem of entropy conservation at the par- 
tonic/hadronic transition with quark recombination. We 
find that the hadronization procedure implemented in the 
qMD model is compatible with the necessary increase of 
entropy with time. 

The charged particle ratio fluctuations D, the charge 
transfer fluctuations D u / (dN c h/dy), the baryon number- 
strangeness correlation coefficient Cbs, the charge- 
strangeness correlation coefficient Cqs were then com- 
puted and found to be compatible with their hadronic 
expectation values in the qMD final state, even though 
in line with the expected QGP results in the early stage. 
The vanishing of the initial QGP correlations and fluc- 
tuations occurs at the partonic/hadronic transition and 
was traced back to the recombination-hadronization pro- 
cess. Our results are supported by the agreement of our 
model with various ratios of susceptibilities calculated in 
lattice QCD. 

The investigation of the variances and covariances di- 
vided by the number of charged particles <jxy = [(XY) — 
(X) (Y)) I (N c h) permits to directly access information on 
diagonal and off-diagonal susceptibilities. We show that 
these quantities exhibit a sharp increase from their QGP 
to their hadronic value at the transition. All the sus- 
ceptibilities present a maximum in the vicinity of the 
hadronization time with the notable exception of the 
electric charge susceptibility ((Q 2 ) - (Q) 2 )/(N ch ). We 
stress that the study of these quantities, only taking 
into account the electric charge Q, baryon number B 



and strangeness S (or another equivalent set of conserved 
charges), allows to treat completely the problem of initial 
QGP correlations and fluctuations for a system of light 
quark flavours u, d and s. This work can be extended to 
heavy quark flavours. 

The transition from the initial QGP to the expected 
hadronic correlation and fluctuation values results from 
various effects. The vanishing QGP signal for charged 
particle ratio fluctuations, charge transfer fluctuations as 
well as for the (Jxy coefficients are driven by the decrease 
of the number of charged particles (N c h) at hadroniza- 
tion. For the ratios of susceptibilities, independent of 
(N c h), the vanishing signal is related to the gathering 
of quarks into color neutral clusters due to the interac- 
tion potential used, the rapidity shift of the quark to 
the hadronic cluster rapidity at hadronization, and the 
decay of resonances. The relative importance of these 
effects differs for different ratios. 

The fact that no similar fluctuation signals was ever ob- 
served in the experimental data until now can thus be un- 
derstood as the result of (recombination-)hadronization. 
It would be interesting to see the experimental result 
of the fluctuation and correlation quantities investigated 
here as they are complementary to the charged particle 
ratio fluctuations, already extensively measured. 
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